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Abstract 



Aims. The mobility of H atoms on the surface of interstellar dust grains at low temperature is still a matter of debate. In dense clouds, 
the hydrogenation of adsorbed species (i.e., CO), as well as the subsequent deuteration of the accreted molecules depend on the 
mobility of H atoms on water ice. Astrochemical models widely assume that H atoms are mobile on the surface of dust grains even if 
controversy still exists. We present here direct experimental evidence of the mobility of H atoms on porous water ice surfaces at 10 K. 
Methods. In a UHV chamber, 2 is deposited on a porous amorphous water ice substrate. Then D atoms are deposited onto the surface 
held at 10 K. Temperature-Programmed Desorption (TPD) is used and desorptions of Ot and D2 are simultaneously monitored. 
Results. We find that the amount of O2 that desorb during the TPD diminishes if we increase the deposition time of D atoms. O2 is 
thus destroyed by D atoms even though these molecules have previously diffused inside the pores of thick water ice. Our results can be 
easily interpreted if D is mobile at 10 K on the water ice surface. A simple rate equation model fits our experimental data and best fit 
curves were obtained for a D atoms diffusion barrier of 22+2 meV. Therefore hydrogenation can take place efficiently on interstellar 
dust grains. These experimental results are in line with most calculations and validate the hypothesis used in several models. 

Key words. Astrochemistry — ISM: Atoms — Dust, extinction — Methods: Laboratory 



1. Introduction 

Among the numerous molecules detected in different astrophys- 
ical envi ronments, a large f raction is composed of hydrogenated 
species (Chang ~et al.ll2.Q07l and references therein). It has been 
established that dust g rains play a major role in the hydro- 
genat ion of compounds ([Herbst & Klemperetll 1973b ISmifh et aT] 
2006). In the dense interstellar medium where most of the hydro- 
genated species have been detected (especially during the star 
formation stage (fTegmark et alJI 19971 : ICazaux & Tielensi r2004; 
Cazaux & Spaans 2004)), grains are covered in a molecular 
mantle mainly composed of water ice. In all the theoretical 
chemical models that describe the hydrogenation of comp ounds 
on grains dTielens & Hagen| [r982; Cup pen & Her bst 2007), mo- 
bility of hydrogen is hypothe sised, though some c alculations 
on amorphous water surfaces ( Smoluchowski 1981) and mod- 
elling of experimental data (Pere ts et al.ll2005l) rule out such a 
mobi l ity. On the other han d, other calculations dBuch & Zhangl 
1991; Masuda et al. 1998) and interpretations of experimental 
data ( Hornekasr et alj2003l) validate the hypothesis that H is mo- 
bile at 10 K, a temperature supposed to be close to that of dust 
grains in interstellar dense clouds where hydrogenation takes 
place. We have already demonstrated that D atoms are mobile at 
10 K (in our experiments tim e scale) on non poro us Amorphous 
Solid Water (np-ASW) ices dAmiaud et alj|2007l) . but the poros- 
ity and the roughness of the water ice films grown at low tem- 
peratures (10 K) may considerably lower the D atom mobility. 

Yet the debate upon the mobility of hydrogen atoms is not 
closed, and i s widely nourished with experim ental interpretation 
controversy dVidali et alJliOOol lDulieul l2006). If one considers 
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only the latest estimated value of 5 1 meV dPerets et al.ll2005l) for 
the diffusion energy barrier of H on porous amorphous solid wa- 
ter (p-ASW) ice, the hopping time between two adjacent adsorp- 
tion sites is about 10 million years at 10 K. In such a context, no 
hydrogenation chemistry can occur on icy mantles of dust grains 
on a reasonable time scale. Hence, observational evidence of hy- 
drogenated and deuterated species is a good argument for an ex- 
perimental investigation of the mobility of H-atoms. 

The aim of this paper is to provide new experimental facts 
that directly address — without the prism of a sophisticated 
model — the question of hydrogen mobility on p-ASW ice at 
10 K. By using O2 as a tracer of D mobility we present a set of 
experiments that are straightforward to interpret if D is mobile 
on the surface of porous amorphous water ice, as it is expected 
from almost all calculations and included in several astrochem- 
ical models. In Section [2] we briefly describe the experimental 
set-up and procedures. In Section|3]we present our experimental 
results and explain them assuming that D is mobile on p-ASW 
ice at 10 K. In Section [4] we describe a simple rate equation 
model that we used to fit our experimental results. In Section 
we discuss other interpretations before concluding. 

2. Experimental Section 

2.1. Experimental Set-up 

The FORmation of MOLecules in the InterStellar Medium 
(FORMOLISM) experimental set-up has been developed with 
the purpose of studying the reaction and interaction of atoms and 
molecules on surfaces simulating the dust grains under interstel- 
lar conditions (relevance of substrate, low density, and very low 
temperature ~ 10 K). FORMOLISM is composed of an ultra- 
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high vacuum chamber with a base pressure of ~1(T 10 mbar, 
of a rotatable quadrupole mass spectrometer (QMS) and of an 
oxygen-free high-conductivity copper sample holder. The sam- 
ple holder is attached to the cold finger of a closed-cycle He 
cryostat and can be cooled down to 8 K. The temperature is 
measured with a calibrated silicon diode clamped on the sample 
holder and controlled by a Lakeshore 334 controller to +0.2 K 
with an accuracy of +1 K in the 8^-00 K range. Reactants are 
introduced into the vacuum chamber via two separated triply dif- 
ferentially pumped beam lines aimed at the surface. The first one 
is used to deposit O2 while the second beam line is used to in- 
troduce the D atoms. This second beam line consists, in its first 
stage, of an air-cooled quartz tube surrounded by a microwave 
cavity for dissociating D2. To cool the D atoms exiting the disso- 
ciation region, the source is terminated with an aluminum nozzle 
connected to a closed cycle He cryostat that can be cooled down 
to ~ 25 K. For the experiments presented here the D beam has a 
temperature of 50 K and a flux of ~1 10 I3 crrT 2 . The dissociation 
rate is constant at ~50%. 



2.2. Growth of the ASW films 

For this work we grow composite ice films of p-ASW on np- 
ASW. The sub-layer of np-ASW has a thickness of -50 ML 
(1 ML = 10 15 cm -2 ) and is grown at 120 K with a micro-channel 
array doser. The only purpose of the compact ice sub-layer 
is to isolate the en suing p-ASW films from the Cu surface 
dEngquist et al.lll995l) . On top of the np-ASW substrate, the p- 
ASW films are grown at 10 K. In the present study we use two 
different thickness of p-ASW, namely 20 ML and 250 ML. The 
20 ML p-ASW film is grown from a background flux formed 
by controlling the H2O partial pressure in the vacuum cham- 
ber. Since ASW films deposited from ambient vapour take a 
long time to form and since a long time is needed to reach the 
base pressure again, the 250 ML p-ASW films are deposited 
with the micro-channel array doser. A previous calibration has 
shown that under these conditions we grow 0.33 ML of ASW per 
second. Although p-ASW films grown by this method exhibit 
a lower degree of porosity, the larger thickness ensures a more 
complex pore structure than that of thin p-ASW ices grown by 
background dosing. After the growth of the composite ice films, 
they are annealed to 70 K to avoid any subsequent morpholog- 
ical changes. Although the porous structure is not completely 
destroyed by the annealing pro cess, this reduces irrev ersibly the 
porosity of the p-ASW films dKimmel et al .112001 alibi) . This is 
verified by looking at the shift in temperature of the D2 desorp- 
tion peak as described in Hornekaer et al. 120051 

2.3. Measurement Procedures 

Our experimental procedures are the following: 

- After the preparation of a stable and pure O2 beam, the p- 
ASW surface held at 10 K (or 25 K in another set of exper- 
iments) is exposed to 0.5 ML of O2 molecules. The flux has 
been previously calibrated using TPD by determining the O2 
exposure t ime required to satur ate the O2 monolayer on np- 
ASW ice (iKimmel et al.ll2001bl) . 

- The p-ASW substrate kept at 10 K is then exposed to vary- 
ing amounts of cold D atoms (50 K) ranging from exposed 
monolayer (EML) to 2 EML. We use EML and not directly 
ML because in the case of light molecules the sticking coef- 
ficient is not necessarily unity. Thus less molecules are ad- 
sorbed on the surface than what the surface has been exposed 
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Figure 1. TPD profiles of D2 and O2 after deposition of 0.5 ML 
of O2 and followed by different deposition times of D (20 sec- 
onds (squares), 40 seconds (triangles) and 90 seconds (circles)) 
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Figure 2. Normalised O2 (plain circles) and D2 (empty squares) 
TPD peak areas following different D exposures. The O2 nor- 
malisation was made with respect to the TPD peak area corre- 
sponding to 0.5 ML of O2 alone and the D2 normalisation was 
made with respect to the TPD peak area corresponding to 1.8 
EML of D. O2 and D both deposited at a surface temperature of 
10 K. The O2 and D both deposited at a surface temperature of 
10 K. The solid line represents the model results (see text). 



to. Moreover in the case of p-ASW ice more than 1 ML of 
D2 can be adsorbed on the surface, and still the coverage can 
be low, since the effective surface ar ea is very large due to 
the 3D conformation of the substrate ( Ami aud et al.l l2006). 

- Finally, the QMS is placed in front of the surface and TPD 
mass spectra of D2 and O2 are recorded simultaneously. All 
TPDs are done with a linear heating rate of 10 K/min starting 
from 10 K up to a temperature above which the O2 molecules 
have completely desorbed (~65 K). 

- The same series of experiments are also done by depositing 
D and D2 followed by O2. 
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3. Results 

3.1. 20 ML p-ASW substrate 

As a first result we find that O2 is only sensitive to D and not to 
D2. When D is deposited on the surface and then exposed, im- 
mediately after that, to O2, we find that the amount of desorbed 
O2 is constant (within the experimental error bar (<10%)) irre- 
spective of the D amounts. This series of experiments was done 
in order to rule out the possibility that D atoms were mobile 
only during the heating phase of the TPD. Indeed these results 
demonstrate that D atoms that adsorb on the ice surface diffuse 
immediately and either react with each other or desorb before O2 
is deposited. The amounts of desorbing O2 are the same regard- 
less of the amounts of D atoms deposited on the surface. Hence, 
in what follows, we restrict ourselves to experiments where O2 
is deposited first and then followed by various D exposure times. 

Figf2]displays the normalised TPD peak areas of O2 and D2, 
as a function of D exposure time (lower axis) or D EML (upper 
axis). FigQ]and Fig|2]show clearly that the quantity of desorb- 
ing O2 decreases with the increase of the exposed dose of D. 
A sudden change in the slope of the monotonic decrease occurs 
after ~ 0.5 EML of D exposure. At first, the drop of O2 that 
desorb from the surface occurs rapidly, then it becomes much 
slower. On the other hand, as expected, the amount of D2 desorb- 
ing from the p-ASW surface increases with the exposure time of 
D (Fig|2|i. Yet, in this case too, we observe that a change in the 
slope of the D2 desorption raise takes place almost exactly after 
~ 0.5 EML of D exposure. Here the increase of D2 TPD peak 
areas is slightly accelerated after the slope change . 

These observations can be explained as follows: as D atoms 
are deposited onto the p-ASW sample, two competing mech- 
anisms occur, namely O2 destruction and D2 formation. If D 
atoms are mobile, they are likely to encounter either another D 
atom to form D2 (D+D formation) or an O2 molecule and react 
with it (O2 destruction). It is significant that in both cases of O2 
and D2 desorptions, the change in slope occurs after the same 
amount of D exposure time (see arrow in Figf2]i. This suggests 
that after ~ 50 s (0.5 EML) of D irradiation, when most O2 has 
been destroyed, a greater number of D atoms are available to 
form D2. The presence of a rapid regime of O2 destruction (80% 
of initial dose destroyed by 0.5 EML of D) and of a slow regime 
(for the remaining 20% of O2 initial dose) can be ascribed to the 
non perfect overlapping of the two beams. It has been checked 
that the two beams have in fact ~ 80% overlap on the surface. 



3.2. 250 ML p-ASW substrate 

We then did the same experiment using a thicker p-ASW ice film 
in order to check whether the destruction of O2 on a 20 ML p- 
ASW ice substrate is actually due to mobile D atoms being able 
to scan the porous surface. This time, on 250 ML p-ASW films, 
we deposited O2 at a higher temperature (25 K), to favour O2 
mobility and to have O2 not only adsorbed on the surface of the 
film but also deeper into the porous structure. The results are 
shown in Fig|3] Plain circles represent the yield of desorbed O2 
after cooling the surface down to 10 K and then irradiated with 
0, 50, 100 and 150 seconds of D exposure at 10 K, while empty 
circles show the case of desorption of O2 initially deposited at 
10 K and then exposed to the same D exposure times. First, we 
notice in these experiments that the O2 destruction rate is much 
slower than in the case of a 20 ML p-ASW substrate. We then 
only observe the first part of the curve in Fig. 2 at a slower rate. 
This can be due to the fact that the roughness is greater in the 
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Figure 3. Comparison between normalised TPD peak areas of 
0.5 ML of 2 deposited at 10 K (blank dots) and 25 K (black 
dots) and then irradiated with 0, 50, 100 and 150 seconds of D 
atoms. The solid line represents the model results (see text). 



case of a 250 ML substrate. The structures (holes, pores, piles...) 
are more pronounced. The total accessible surface is also higher 
and subsequently the average density of O2 is reduced. In this 
case D2 formation is more favoured. Secondly, we find that if O2 
is deposited at 25 K, the destruction rate is even lower. This can 
be due to the fact that at 25 K O2 is more mobile on the surface 
of the ASW ice than at 10 K. Therefore O2 molecules are spread 
over all the accessible porous surface. O2 density is again lower 
and D2 formation is favoured. We also notice the disappearance 
of the elbow that was present in Fig. 2 at 0.5 EML of D. It is 
due to the fact that we should have exposed the surface to longer 
times of D atoms in order to reach 80% of destruction of the 
adsorbed O2 to observe it. Finally we have checked that all the 
O2 deposited at 25 K disappears, after very large exposure times 
(600 seconds), even on the 250 ML p-ASW film. 

4. Model and data fit 



A very simple model based upon the formalism of Kat z et ail 
( 1999) was used to fit the experimental results already obtained 
and shown in Fig. 2. If F tol is the initial beam flux of D2, after the 
dissociation F, of is divided into a flux of D, F/j, and a flux of D2, 
Foi- Both fluences on the surface are determined knowing the 
sticking probabilities of D2 and of D, Sdi and So respectively, 
and the dissociation rate rat a beam temperature of 50 K. The 
fluxes are given by: 

F/j = 2rP D F, , and F D2 = (1 - r)P D2 F fof 

Taking into account the porosity of the ice surface, the surface 
density of D, D2 and O2, and depending on the diffusion barrier 
Erf/// of D atoms on the surface, this model simulates the number 
of D + D reactions as well as the number of D + O2 reactions. As 
Fig. 2 and 3 clearly show, our simple model is able to reproduce 
the experimental results under the same conditions, namely, dis- 
sociation efficiency of the D2 beam, sticking probabilities of D 
and D2 (measured at a given beam temperature, Matar et al. in 
preparation), and porosity of the ice. A best fit of the experimen- 
tal data has given a diffusion barrier of 22+2 meV. 
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5. Discussion 

In the previous section we have proposed a direct interpretation 
of our experimental data assuming that D is mobile at 10 K. But 
other hypothesis should be also discussed in the light of these 
new experimental data, in particular the possible thermally in- 
duced mobility of the atoms. This process implies that atoms are 
not mobile at 10 K and that they require higher temperature to 
move and react. If atoms are not mobile at 10 K, they should stay 
in the vicinity of the external surface, especially under our exper- 
imental conditions where D atoms have a kinetic temperature of 
50 K (in the beam) which reduces considerably the mobility of D 
atoms resulting from the hot atom (Harris-Kasemo) mechanism. 

The total number of D atoms sent to the surface can be quite 
high (2 10 cm -2 , i.e. 2 ML) which corresponds to a value of 
full coverage of the external surface. Consequently, if atoms are 
not mobile, the surface density becomes high, and two mech- 
anisms invoked in the literature should become non negligible. 
One is the direct Eley-Ridel formation mechanism for D2 (an 
adsorbed atom forms a molecule with an atom incoming from 
the gas phase, a fraction of the formed molecules is then directly 
released in the g as phase), the second is the rejection mechanism 
dPerets et a l. 2005) (adsorbed atoms prohibit atoms from the gas 
phase to stick on the surface). In both cases a fraction of D or D2 
evaporates during the exposure and therefore does not desorb 
later on during the heating ramp. If so, D2 should not be pro- 
portional to the exposure, which is not the case in Fig. 2 where 
we observe proportionality. We can thus conclude that either the 
surface density does not increase to reach a value of the cover- 
age close to unity, either the two cited mechanisms (Eley-Ridel 
and rejection mechanism) are not efficient. It is really unlikely 
that the Eley-Ridel mechanism does not occur, therefore density 
would not increase, implying that D atoms enter the pores. 

Calculations dBuch & Czerminskllll991l) show that the ad- 
sorption energy of H on a water cluster and the height of the 
energy barrier between two adjacent sites are both described 
by large distributions. This means that the mobility is site de- 
pendent. Some deep adsorption sites can bind atoms more than 
others and release them less easily in the gas phase. Thus it is 
possible that a fraction of the atoms are not mobile at 10 K, 
on the time scale of the experiment, and require thermal acti- 
vation to react. Our experiments just show that globally the D 
atoms are able to penetrate in the porous structure of the ice. 
In fact, we find a non-infinite D atom mobility E^yy = 22+2 
meV, and even if the model is based upon pure thermal hopping, 
this diffusion barrier physically means that the mean residence 
time before hopping is ~12 ms. Our results do not show if this 
mobility is dominated by thermal hopping or tunnelling. More 
experiments should be done at lower surface temperature (be- 
low 8 K) to fully understand how this mobility occurs. Another 
aspect of these experiments concerns the efficiency of the O2 + 
D reaction in comparison with the D + D reacti on. We can see 
that O2 + D — » O2D is a very efficient reaction (Mivauchi et al.l 
2008; Ioppolo^etal. , 2008), because the proportion of D required 
to destroy O2 is close to 1, at the accuracy of our fluxes esti- 
mation (~50% ), because O2+D is believed to have a low en- 
ergy barrier dWalch et al. II 1988b . In our experiments the fact that 
the O2+D reaction seems to be favoured in comparison with the 
D+D reaction when the exposition is low does not signify that 
the D+D reaction is less efficient because of an activation barrier. 
Indeed, our model assumes that both D+D and O2+D reactions 
have 100% efficiency. Once again, if atoms are mobile, the den- 
sity of D atoms stays low because they are consumed by the O2 
population. It is only when the density of D atoms is equal to 



that of O2 atoms that the production of D2 can be effective. This 
consideration allows us to conclude that at any moment in our 
experiment, the density of D atoms is low, initially because they 
are consumed by O2 and subsequently because they form D2. 

At the end of the experiment, when the ice is sublimated we 
detect a high prop ortion of HDO and not D2O because of the 
isotopic exchange (Smi th et"aT1ll997l) . Therefore O2 is one of 
the possible precursor for water formation. Finally these experi- 
ments show that no O2 should survive on grains in dark interstel- 
lar clouds. Apart from specific detection limitations, it can be a 
fundamental reason of the non detection of O2 on icy grains in 
dark clouds. 

6. Conclusions 

We have reported the results of experiments on the mobility of 
deuterium at 10 K on porous ASW ice surfaces under interstel- 
lar conditions, using a temperature programmed desorption tech- 
nique. Beams of O2 and D were irradiated on the surface of an 
p-ASW ice film and the mobility of D atoms at 10 K was inves- 
tigated via their property of reacting with the O2 molecules. By 
using a simple rate equation model we found that the D atoms 
diffusion barrier at 10 K is 22+2 meV. Our experimental results 
demonstrate that D atoms are mobile on ASW ice surfaces at 
10 K, and validate the chemical models that hypothesized the 
mobility of H atoms at 10 K. 
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